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ABSTRACT 

We describe a novel technique for probing the statistical properties of cosmic magnetic fields 
based on radio polarimetry data. Second-order magnetic field statistics like the power spec- 
trum cannot always distinguish between magnetic fields with essentially different spatial 
structure. Synchrotron polarimetry naturally allows certain 4th-order magnetic field statistics 
to be inferred from observational data, which lifts this degeneracy and can thereby help us 
gain a better picture of the structure of the cosmic fields and test theoretical scenarios describ- 
ing magnetic turbulence. In this work we show that a 4th-order correlator of specific physical 
interest, the tension-force spectrum, can be recovered from the polarized synchrotron emis- 
sion data. We develop an estimator for this quantity based on polarized-emission observations 
in the Faraday-rotation-free frequency regime. We consider two cases: a statistically isotropic 
field distribution, and a statistically isotropic field superimposed on a weak mean field. In both 
cases the tension force power spectrum is measurable; in the latter case, the magnetic power 
spectrum may also be obtainable. The method is exact in the idealized case of a homogeneous 
relativistic-electron distribution that has a power-law energy spectrum with a spectral index 
of p = 3, and assumes statistical isotropy of the turbulent field. We carry out numerical tests 
of our method using synthetic polarized-emission data generated from numerically simulated 
magnetic fields. We show that the method is valid, that it is not prohibitively sensitive to the 
value of the electron spectral index p, and that the observed tension-force spectrum allows 
one to distinguish between, e.g., a randomly tangled magnetic field (a default assumption in 
many studies) and a field organized in folded flux sheets or filaments. 

Key words: galaxies: clusters: general; intergalactic medium; ISM: magnetic fields; magnetic 
fields; methods: data analysis; radio continuum: general; turbulence 



1 INTRODUCTION 

Magnetized plasma is present almost everywhere in the observable 
Universe, from stars and accretion disks to the interstellar and the 
intracluster medium (respectively ISM and ICM). A large fraction 
of this magnetized plasma is in a turbulent state. Understanding 
the origin of the cosmic magnetic fields and their evolution to- 
wards their observed state embedded in magnetized plasma turbu- 
lence, apa rt from being a tantalizing int ellect ual challenge in its 



own ri ght dBrandenburg & Subramania n 2005; Subramanian et al 
ISchekochihin & Cowlevl 120071 120061 : ISchekochihin et al 



2006 



2009), is also crucial in the construction of theories of large-scale 
dynamics and transport in many astrophysical systems. For exam- 
ple, magnetic fields are expected to be dynamically important in 
determining the angu l ar momentum transport in accretion discs 
( tPringle & Reesll 19721 : IShakura & Svunvaev|[l973l) . to control star 
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formation and the general structure of the ISM (where magnetic 
fields prevent molec ular clouds from collapsing and suppress frag- 
mentation; see, e.g.. |Price et al .1120081 and references therein), and 
to play an important role in galaxy discs as well as galaxy clusters, 
where they infl uence the viscosity and thermal conductivity of the 
ISM and ICM dChandran & Cowlevl [l 9981 : iNaravan & Medvedevl 
l200ll:lMarkevitch et alj|2003[) and the propaga tion of cosmic rays 
(e.g JStrong et alj|2007l: lYan & Lazarianll2008l) . Theoretical mod- 
els of all these phenomena require some assumptions to be made 
about the spatial structure of the tangled magnetic fields permeat- 
ing the constituent turbulent plasmas. However, as theory of mag- 
netized plasma turbulence is in its infancy as a theoretical sub- 
ject, there is no consensus about what this spatial structure is. In 
order to make progress both in understanding the turbulence and 
in modeling its effect on large-scale dynamics and transport, it is 
clearly desirable to be able to extract statisti cal information about 
the field structure from observational data i lEnfilin& Vo"g3l2006l : 
lEnfflin et al.ll2006T) . 
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Figure 1. Left panel: MHD-dynamo-generated magnetic field {Schekochihin et al. 2004j) (saturated state of their Run S4). Right panel: A synthetic divergence- 
free Gaussian random field realization with identical power spectrum. These are cross-sections of the field strength |B| (dark represents stronger field, white 
weaker field). The magnetic-field and tension-force power spectra are shown in Fig. [2] 



Diffuse synchrotron emission is observed throughout the 
ISM and the ICM, as we l l as in the lobes of radio galaxies 



1994; iReichet al 


|200l|; 


Beck et alj |2002|; IWolleben et alj 12006]; 


Haverkorn et alj 


2006 


\ 


Reicbl 120061; Clarke & EnMinT 120061; 


Schnitzeler et alj 


2007 


; lLains et al. 20081) 


The fact that syn- 



chrotron emission is readily observable and is a good tracer of the 
magnetic-field strength and orientation makes it a key source of 
information that can serve as a reality check for theories of mag- 
netized plasma turbulence and magnetogenesis (origin of the mag- 
netic fields). 

In this work, we will be focusing on how the synchrotron- 
emission data can be used to characterize the structure of the tan- 
gled magnetic fields permeating the ISM and the ICM. In this con- 
text we refer to previous studies which sought to recover statis- 
tical information about the structure of of th ese fields in the ICM 
from the Faraday rotat i on me a sure (RM) data (EnBlin & Voga2003l; 

ovoni et al] 120061 ; iGuidetti et alj 



trom the Faraday rotation measure 
Vogt & EnBlir] l200j 120051 ; |Gc 



20081) . as well as studies of the ISM dHaverkorn et alj|2006ll2008l) . 
als o based on the RM data, and the work of lSpanglerl ( 1 1983. 1 19831) 



and lEilekl ( fl989al lrjh based on polarized synchrotron emission data. 
In formal terms, all of these papers are concerned with at most 
second-order statistics, namely the magnetic-field power spectrum, 
or the two-point correlation function of the magnetic field. Our 
work complements those previous efforts by drawing on the fact 
that polarized-emission data carries information about 4th-order 
statistics of the magnetic field. In particular, we present a practical 
method for obtaining the tension-force power spectrum. As will be 
shown in greater detail in the following, this quantity contains sta- 
tistical information about the spatial structure of the tangled mag- 
netic fields that is missing in the second-order statistics and, most 
importantly, is actually observable with radio telescopes mapping 
polarized synchrotron emission. 

The plan of this paper is as follows. In Sec. [2] we explain why 
the tension-force power spectrum is an interesting quantity to mea- 



sure and how it allows one to diagnose the magnetic-field structure. 
In Sec. [3] we explain the assumptions we make about the magnetic 
field (Sec. 13.1 1 and Sec. |3.2| > and the observational data (Sec. 13.31 
see also Appendix [Aj and propose a method of reconstructing the 
tension-force power spectrum from the Stokes maps (Sec. |3.4| and 
Sec. l3.5t , We then generalize our method slightly for the case when 
a weak mean field is present and show that in this case the power 
spectrum of the magnetic field itself may be obtainable from the 
Stokes correlators (Sec. l3.6t . Most detailed analytical calculations 
required in this section are exiled to Appendix [B] In Sec. [4] we 
demonstrate the validity of our method by testing it on synthetic 
observational data generated from numerical simulations. A brief 
summary and conclusion is given in Sec. [5] 



2 MOTIVATION 

Turbulent plasmas exhibit in general very complex magnetic struc- 
tures (see Fig.[T] left panel), which are best characterized by statis- 
tical means. The most widely used quantity for this purpose is the 
power spectrum 



M(k) = 4nk 2 {\B(k)\ 2 ), 



(1) 



where B(k) is the Fourier transform of the magnetic field (see 
Sec. 13. It . The angle-bracket averaging includes averaging over all 
directions of k, so the power spectrum measures the amount of 
magnetic energy per wavenumber shell \k\ = k. It is related via 
the Fourier transform to the second-order two-point correlation 
function (or structure function) of the magnetic field. It is an at- 
tractive quantity to measure because phenomenological theories of 
turbulence typically produce predictions for characteristic field in- 
crements between two points separated by a distance in the form 

of power laws with respect to that distance dKolmogorovl 
1 1 I 1 I 1 1 I -i I ^^^^^^^^^^^ 

Iroshnikov 1963; Kraichnan 1965; Gold reich & Sridharl 

Boldvrevl 120061 ; ISchekochihin et alj 120091) — and such predictions 



194ll; 



1995; 
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Figure 2. Left panel: Magnetic-field power spectra for the fields shown in Fig.fT] Right panel: Tension-force power spectra for the same fields. 



are most obviously tested by measuring the spectral index (or the 
scaling exponent of the structure function). However, knowing the 
spectrum is not enough and can, in fact, be very misleading, for 
reasons having to do both with the physics of magnetic turbulence 
and with formal aspects of describing it quantitatively — the basic 
point, which is discussed in much detail below, being that spectra 
do not contain any information about the geometrical structure of 
the magnetic-field lines. 

All scaling predictions for magnetized plasma turbulence pro- 
posed so far are, implicitly or explicitly, based on the assump- 
tion that magnetic fluctuations at sufficiently small scales will 
look like small Alfvenic p erturbations of a larger-scale "mean" 
field (this is known as the iKraichnanl [l965l hypothesis). Numer- 
ical simulations of MHD turbulence carried out without impos- 
ing such a mean field do n ot appear to support this hypothe- 
sis dSchekochihin et alj|2004h . although the currently achievable 
resolution is not sufficient to state this beyond reasonable doubt 
and t he results are to some extent open to alternat ive interpreta- 
tions faaugen et aT]|2004 1 Subramanian et al. I l2006|) . What seems 
to be clear is that the magnetic field has a tendency to or- 
ganize itself in long filamentary struct ures ("folds") with field- 
direction reversals on very small sc ales dSchekochihin et al.ll2004l : 
iBrandenburg & Subramanianll2005t) . Filamentary magnetic struc- 
tures are, indeed, observed in galaxy clusters dEilek & Owen! 
|2002| ; IClarke & Enfiiinl l200ej) . although the field reversal scale 
does not appear to be nearly as small as implied by MHD 
turbulence simulations — a theoretical puzzle solving which will 
probably require bringing in kinetic physics (see discussion in 
ISchekochihin & Cowley||2006l) . 

It is clear that both the current and future theoretical debates 
on the structure of magnetic turbulence would benefit greatly from 
being constrained observationally in a rigorous way. For the rea- 
sons explained above, in order to do this, we must be able to diag- 
nose nontrivial spatial structure, which cannot be done by looking 
at the magnetic power spectrum alone. Let us explain this in more 
detail. 

Consider a divergence-free, helicity-free, statistically homo- 
geneous and isotropic field as a minimal model for the fluctuating 
component of the magnetic field in galaxies and clusters. If this 
field also obeyed Gaussian statistics exactly or, at least, approxi- 
mately, its power spectrum would be sufficient to completely de- 
scribe its statistical properties because all higher-order multi-point 
statistics could be expressed in terms of the second-order two-point 



correlators and, th erefore, the power sp ectru m. As s uming such 
Gaussian statistics. Ispanglerl d 19821 1 19831) and lEilekl dl989ah pro- 
posed to calculate the magnetic power spectrum using the observed 
total and polarized synchrotron radiation intensity, quantified by 
the Stokes parameters /, Q and U (see Appendix [A}. Computing 
two-point correlation functions of the Stokes parameters, hence- 
forth referred to as Stokes correlators, one essentially obtains two- 
point, 4th-order correlation functions of the magnetic field in the 
plane perpendicular to the line of sight (Sec. |3.4t . If the statistics 
are Gaussian or if Gaussianity is adopted as a closure assumption, 
the 4th-order correlators can be split into second-order correlators, 
so the power spectrum follows. 

The problem with this approach to magnetic turbulence is that 
the Gaussian closure essentially assumes a structureless random- 
phased magnetic field, which then is, indeed, fully characterized 
by its power spectrum. It is evident in Eq. l[T) that all phase in- 
formation, which could tell us about the field structure, is lost 
in the power spectrum. As we explained above, both numerical 
and observational evidence (and, indeed, intuitive reas oning; see 
ISchekochihin et alj2004llSchekochihin & Cowley|2007t) show that 
magnetic fields do have structure and are very far from being a 
collection of Gaussian random-phased waves. Their spectra tell us 
little about this structure. This rather simple point is illustrated in 
Fig. [T] the right panel depicts an instantaneous cross section of a 
3D magne tic field obtained in a typic al MHD dynamo simulation 
taken from lSchekochihin et alj d2004l) . while the left panel shows a 
synthetically generated divergence-free Gaussian random field with 
exactly the same power spectrum (shown in Fig. [2] left panel). The 
folded structure discussed above is manifest in the simulated field 
but absent in the Gaussian one: in the former case, the field typi- 
cally varies across itself on a much shorter scale than along itself 
and the regions of strongest bending are well localized, whereas in 
the latter case, the field is uniformly tangled and has similar varia- 
tion along and across itself. 

So how can one differentiate between such different fields in a 
systematic and quantitative way (i.e., ot her than by simply looking 
at vis ualizations)? As was pointed out bv lSchekochihin et al.ld2002l 
this can be done by looking at the statistics of the tension 

force 



F = 



B VB 

4tt ' 



(2) 



As a formal diagnostic, the tension force is a measure not just of 
the field strength but also of the gradient of the field along itself, 
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thus it is strong if a field line is curved, and weak if the field line 
is mostly straight. The tension-force field associated with a folded 
magnetic field (strong, straight direction-alternating fields in the 
"folds", weak curved fields in the "bends") will obviously be very 
different from the one associated with a random Gaussian field. As 
shown in Fig.[2](right panel), their power spectra 



T(k) = Ank 2 (\F(k)\ 2 ) 



(3) 



do, indeed, turn out to be very different: flat for the folded field, 
peaked at the smallest scales for the Gaussian field. Why a flat 
tension-force spectrum is expected for a folded field is discussed in 
Schekochihin et al. I I2004 their § 3.2.2, where numerical measure- 
ments of the tension-force statistics can also be found. In contrast, 
for the Gaussian field, one obviously gets T(k) cc k 2 M 2 , hence the 
peak at the small scales. 

In physical terms, the tension force is one of the two compo- 
nents of the Lorentz force 



1 



.B 2 



-JxB = -V— + 

c 8tt 



B VB 



An 



(A) 



where the first term on the right-hand side is the magnetic pressure 
force, and the second term is the magnetic tension force, as defined 
in Eq. {2j. In subsonic turbulence, the tension force essentially de- 
termines the dynamical back reaction of the magnetic field on the 
plasma motions because regions with higher magnetic pressure can 
be expected to have correspondingly weaker thermal pressure, so 
that the magnetic pressure forces are mostly balanced by oppositely 
directed thermal pressure forces. 

Thus, measuring tension-force power spectra not only permits 
one to discriminate quantitatively between different magnetic tur- 
bulence scenarios but also provides a detailed insight into the MHD 
physics occurring in space, because it quantifies the properties of 
the dynamically relevant force in the magnetic turbulence. It is per- 
haps worth stressing this last point. In principle, many 4th-order 
statistical quantities that one might construct out of the Stokes cor- 
relators should be able to discern between different magnetic-field 
structures, but the tension-force power spectrum also has a clear 
physical meaning. 

It is a stroke of luck that not only the tension-force power spec- 
trum is the diagnostic that we would ideally like to know from the 
theoretical point of view, but it turns out that, under mild simpli- 
fying assumptions, it can be fully recovered from the statistical in- 
formation contained in the Stokes correlators and, therefore, it is 
observable! This will be demonstrated in detail in the following 
sections. Such an outcome is not automatic: other potentially in- 
teresting statistical quantities such as the magnetic-energy power 
spectrum or the magnetic pressure-force statistics are not so di- 
rectly imprinted into the Stokes correlators and require further as- 
sumptions in order to be extractable from the same data. 



3 METHOD 

In this section, we outline a formal theoretical framework for con- 
verting polarized-emission observables into the physically interest- 
ing statistical characteristics of the magnetic field under a number 
of simplifying assumptions. 



3.1 Magnetic Field 

Let us assume some volume V of interstellar or intracluster plasma 
to be filled with a magnetic field B(x) and a magnetized relativistic 



electron population giving rise to the synchrotron emission we ob- 
serve (Fig.[3}. We use a Cartesian coordinate system (x, y, z), where 
z is the line of sight. The volume under consideration is assumed 
to have depth L in this direction. The magnetic field can be decom- 
posed into two parts: 



B = B + b, 



(5) 



where B = (B) is the regular (mean) field throughout the volume 
under consideration and b is the fluctuating ("turbulent") field. The 
former is assumed to be known and the latter is what we aim to 
study. We will work out its various correlation functions and their 
relationship to observable quantities — this can be done both in po- 
sition space and in Fourier space in largely analogous ways. The 
Fourier transform of the field is defined according to 



b(k) 



- f d 3 xe- ikx b(x), b(x) = Y J eik ' x kk), 



(6) 



where x = (x, y, z). In what follows we will drop the hats on the 
Fourier transformed quantities. Note that discrete and continuous 
wave- vector spaces are related via a simple mnemonic: 



y «, J_ f d i k . 



(7) 



3.2 Assumptions: Homogeneity and Isotropy 

We will make two key assumption about the fluctuating magnetic 
field: statistical homogeneity and isotropy. The first of these is not a 
serious restriction of generality as, essentially, we would like to cal- 
culate statistical information based on data from subvolumes within 
which system-size spatial variation of the bulk properties of the as- 
trophysical plasma under consideration can be ignored. The second 
assumption, the isotropy, is more problematic because of the known 
property of magnetized turbulence to be strongly anisotropic with 
respect to the direction of the mean field, provided the mean field 
is dynamically strong (see discussion and exhaustive refere nce lists 
in lSchekochihin & Cowley|2007l : ISchekochihin et al.l2009h . It will, 
therefore, only be sensible to apply our method to astrophysi- 
cal situations where the mean field is either absent or weak, i.e., 
B <k (\b\ 2 ). This should be a very good approximation for the 
ICM and may also be reasonable in p arts of the ISM (e.g., in the 
spiral arms; see lHaverkorn et al J 20061 [20081) . 

In what follows, we will first consider the case of B = and 
then provide a generalization of our results to the case of a weak 
mean field (Sec. 13. 6t . In both cases, we will first show how far 
one can get without the isotropy assumption and then find that only 
assuming isotropy are we able to calculate the tension-force power 
spectrum. It will also turn out that, in the case of a non-zero weak 
mean field, additional information can be gleaned from polarized- 
emission data, including the power spectrum of the fluctuating field 
(normally not available without the Gaussian closure, as discussed 
in Sec.0- 

Physically, we might argue that a weak mean field does not 
modify the turbulent dynamics and, therefore, does not break the 
statistical isotropy of the small-scale turbulent field. Obviously, if 
the bulk of the magnetic energy turns out to reside above or at 
some characteristic scale l B , the statistically isotropic fluctuating 
field at that scale will look like a (strong) mean field to fluctuations 
at scales smaller than l B and assuming isotropy of those fluctuations 
will almost certainly be wrong. Thus, our method can only be ex- 
pected to handle successfully magnetic fluctuations at scales larger 
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than 1 B . This, however, is sufficient to make the outcome interest- 
ing because the key question in the theoretical discussions about 
the nature of the cosmic magnetic turbulence referred to in Sec. [2] 
is precisely what determines 1 B (is it the reversal scale of the folded 
fields? what is that scale?) and how diagnosing the spatial structure 
of the field at scales above I B might help us answer this question. 

Note that a field organized in folds or filaments, as in Fig. Q] 
(left panel), is statistically isotropic because, while the folds extend 
over long distances, their orientation is random. 

3.3 Observables: Stokers Parameters 

Our direct observable is the partially linearly polarized synchrotron 
emission of the relativistic electrons gyrating in the magnetic field. 
This emission is measured by radio telescopes in projection onto 
the sky in terms of the Stokes parameters /, Q and U. Let us briefly 
recapitulate the relevant physics. 

We assume a relativistic electron population that is spatially 
homogeneous, has an isotropic pitch-angle distribution, and a 
power-law energy distribution: 



N(y)dy = Cy~ p dy, 



(8) 



where y is the Lorentz factor, N the number of electrons per y per 
unit volume and C is the normalization factor proportional to the 
electron density. The observed emission wi ll then be partially lin- 
early polarized dRvbicki & Lightmarj[l979f) and, therefore, at any 
given observed (radio) frequency v, it is fully characterized by the 
Stokes parameters /, Q and U as functions of the sky coordinate 
x ± = (x, y) (the spatial coordinate in the plane perpendicular to the 
line of sight). This is detailed in Appendix lAl 

As explained in Appendix lAl a measurement of the absolute 
values of /, Q and U (and, therefore, of the magnetic field and 
its tension) requires knowledge of the relativistic-electron energy 
density, which enters via the factor C in Eq. {S). In our Galaxy, 
it can be inferred directly from its value measured at Earth, or 
via independent messengers such as inverse Compton emission of 
starlight and cosmic-microwave-background photons. The absolute 
rms value of the magnetic field can then be inferred either via the 
Stokes parameters or via the popular assumption of equipartition 
between the average magnetic en ergy density and the r elativistic- 
electron energy density (see, e.g. iBeck & Krause|[2005t and refer- 
ences therein). Note that it is essential for the method developed 
below that the relativistic electrons can be assumed to have a homo- 
geneous distribution (C = const), and do not follow the magnetic- 
field-strength enhancements on small scales. Otherwise our method 
would provide incorrect results, since such a local coupling of rela- 
tivistic electrons and magnetic fields is not incorporated, and would 
destroy the assumed relation between the spatial variation of the 
synchrotron observables and the magnetic fields. However, since 
the relativistic electrons are very diffusive along and even perpen- 
dicular to the magnetic fields, a roughly homogeneous distribution 
can safely assumed in most relevant environments. 

We further assume p = 3 in Eq. {8} (corresponding to 
the frequency distribution oc v -1 ). This is a convenient choice 
because then all Stokes parameters are quadratic in the mag- 
netic field, which means that their two-point correlation func- 
tions will give us 4th-order statistics. We stress that this power 
law, although ex pected by theoretical shock acceleration mod- 
els torurvlll983h. is. of course, a simplification of reality (see, 

yjl2 



e.g. IStrong et alj|2007l) . However, it is usually a sufficiently good 
approximation over fairly wide frequency ranges for many syn- 
chrotron sources. Thus, p = 3 is reasonably close to the values ob- 



N 
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Y 
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Figure 3. Magnetic field and the observables: artist's impression. 



serve d for our own Galaxy l lReich &~R eich 1988]; iTatevama et al.l 
Il986h . the v alues obtained by CMB fore g round subtraction tech- 
niques (e.g. iTegmark & Efstathioul 1 19961: Ide Oliveira-Costa et al.l 
l2008tlDunklev et alj2009l ; lBottino et alj2008l). and found in extra- 
galactic observations of radio- galaxies jBeck et alJll996l) . While 
the theoretical developments that follow do depend on taking p = 3, 
the numerical tests of the resulting method reported in Sec. l4.2l will 
show that it is not essential that p = 3 be satisfied particularly pre- 
cisely. Deviations from p = 3 can be addressed analytically in a 
more quantitative way by a Taylor expansion around p = 3, which 
we leave for further work. 

Finally, we assume the observed volume to be optically thin 
and its Faraday depth to be negligible at the observation frequency 
v. At high frequencies, both conditions tend to be satisfied, Fara- 
day rotation being in most cases the greater constraint. For exam- 
ple, in our Galaxy, Faraday rotation is a relevant phenomenon at 
frequencies below a few GHz, while the medium remains mostly 
optically thin down to frequencies of a few hundred MHz, w here 
free-free absorption starts being relevant (see lSun et al .1120081 . and 
references therein). In cases where Faraday rotation is present in the 
frequency range of the data, we assume that a Faraday de-rotation 
has been applied. Even in the case of source-intrinsic Faraday ro- 
tation, this can still be achieved us ing Faraday tomography tech- 
niques terentiens & de Bruvnll200l) . 

Under these conditions the Stokes parameters can be written 
as the following line-of-sight integrals 



x x ) = j f dz[B;(x) + B 2 y (x)], 

*-> Jo 

Q(xJ = ij^ dz[B 2 x (x)-B 2 y (x)], 

U(x ± ) = j f 
L Jo 



dz2BAx)B v (x), 



(9) 



where L is the depth of the emission region. The dimensional pref- 
actors converting the magnetic-field strength to radio emissivity 
have been suppressed (see Appendix lAl. 
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3.4 From Stokes Correlators to Magnetic-Field Statistics 

Thus, observed polarized emission provides us with three scalar 
fields related quadratically to the magnetic field projected onto the 
plane perpendicular to the line of sight. We can construct 6 two- 
point correlators of these fields, which we will refer to as the Stokes 
correlators: 

SxrOx) = (X(x x ) Y(x ± + r x )>, (10) 

where X,Ye {I,Q,U} and (...) denote a statistical average per- 
formed over the observational maps, which usually means volume 
averaging with respect to the sky coordinate x ± . 

Are the Stokes correlators sufficient to reconstruct the statis- 
tics of the magnetic field? 

In formal terms, the statistical properties of a stochastic field 
are fully described by its rc-point distribution function, or, equiva- 
lently, by the full set of its ra-point, m-th order correlation tensors. In 
practice, this is too much information, most of it is not observable 
in any realistic situation, and in any event, only a few of these corre- 
lators can be interpreted in simple physical terms and are, therefore, 
useful for a qualitative understanding of the field structure. As the 
Stokes correlators are 4th order in the magnetic field and measure 
its correlations between two points in space, it is the two-point, 
4th-order correlation tensor that will be relevant to this discussion: 

C,j, m „{r) = (Bi(x)Bj(x)B m (x + r)B„(x + r)> 

= (H^xWUx + r)), (11) 

where, for notational convenience, we have introduced the field ten- 
sor Hjj = BjBj. The angle brackets denote statistical average, un- 
derstood ideally as an ensemble (or time) average and in practice, 
if we are dealing with one observed realization of the field, as the 
volume average: (...) = (1/V) f d 3 x(. . .). Implicitly, performing a 
volume average relies on the assumption of statistical homogene- 
ity (Sec. |3.2| l, i.e., independence of the statistical properties of the 
field of the reference point x where they are calculated. In terms of 
Fourier-space quantities, we have 

C,;,„„,(r) = ^ e ikr C Ujm ,(k\ C iim (k) = (Hjj(k)H m (k)), (12) 

where the Fourier transforms of all quantities are defined similarly 
to Eq. 

In general, the tensor Cjp„„ depends on very many inde- 
pendent scalar functions, so the 6 available Stokes correlators 
[Eq. d lOb l cannot provide all the required information necessary to 
recover the magnetic-field statistics. Indeed, let us write the Stokes 
correlators in terms of the correlation tensor C, p, m . It is particularly 
easy to do this in Fourier space because the line-of-sight integration 
in Eq. {9} amounts simply to picking the k- = component of the 
field: 

/(ft,.) = + Hyy(k ± ), 

Q(k ± ) = H tx (k ± )- H yy (k ± ), 

U(k ± ) = 2H xy (k ± ), (13) 

where k ± = (k x ,k y ,0). Therefore, the Fourier transforms of the 
Stokes correlators [Eq. i ll OH are 

E 7 /(*x) = C xx , xx (k ± ) + C xx ,yy(k ± ) + C XX j, y (k ± ) + C yy , yy (k ± ), 

^QQ(k±) = C xx , xx (k ± ) - C xx ,yy(k x ) - (":, ,,(*. ) + C,,. yy (A ± ), 

Zuu(k±) = 4C„,.„(A X ), 

2/e( fe ±) = C XXrXX (k ± ) - (',,. (A.) + C xxyy (k ± ) - C yytyy (kJ, 

Ziu(k x ) = 2|(".. (A ) - (".. ..(A. )|. 



Z QU (k ± ) = 2|C, (A i C . (A )j. (14) 
This immediately implies that 

C xx , xx (k±) = 4 P;/(A ± ) + X QQ (k ± ) + 2Re E Je (*J] , 

C yy , yy (k ± ) = - p ;/ (A ± ) + l QQ (k ± ) - 2Re £ Je (* ± )] , 

C xx , yy (k ± ) = - [Z,/(k ± ) - £ fie (Aa) - i'2Im S/ e (fc x )] , 

C.vv,vy(Aj_) = —T.uu(k±), 

C <k) = -[S w (fc x ) + S e[r (fcJ], 

C , >A ) = ip/^JfeJ-Sfl^JfcJ]. (15) 

These are the only components of the correlation tensor Eq. J 1 2b 
that are observable directly and it is only their dependence on the 
wave vector perpendicular to the line of sight that can be probed. 
No correlators that involve the projection of the field on the line of 
sight (B-) can be known. 

The number of independent scalar functions that determine 
Cijj„„ is reduced and becomes closer to the number of observables 
if we make some symmetry assumptions and, in particular, isotropy 
(Sec. |3.2t . Under this assumption, it turns out that we only need to 
know 7 independent scalar functions of k = |A| to reconstruct C, JiM „ 
fully (see Appendix IB 1.3\ . It also turns out that, if no mean field 
is present (B = 0), only 4 of the Stokes correlators of an isotropic 
field contain independent information: S//, two of ~Lqq, Zw, ^qu 
and one of ~L IV , For example, if we keep ~L UV and 

S/q, the other two Stokes correlators are 

2/t/ = S/2 tan 2ip, 

T,qu = - (Z.qq - Lf/f/)tan4y>, (16) 

where ip is the angle between A x and the x axis of the 
frame in which the Stokes parameters are measured, i.e., A ± = 
A x (cos^>, sin<^,0). These relations are useful in constructing well 
behaved expressions for the observables (see Sec. 13.51 and Ap- 
pendix |B2.2| (. They could also be useful in practical situations when 
the Stokes maps might not be perfect, so one might have more (or 
higher-quality) data on some Stokes correlators than on others. 

We see that, even with isotropy, we do not have enough ob- 
servables to measure the general 4th-order statistics of the magnetic 
field (7 independent scalar functions needed, 4 available). How- 
ever, the information carried by the Stokes correlators does suffice 
to reconstruct some of the correlation functions of the field. How 
to determine whether any particular 4th-order correlator is observ- 
able is explained in Appendix IB2.2I In a stroke of luck, we find 
that we can reconstruct the tension-force power spectrum, which is 
a physically interesting quantity because it diagnoses the geomet- 
rical structure of the magnetic field and its dynamical action on the 
plasma motions (Sec.O. Although it follows from the general pro- 
cedure given in Appendix |B2.2| (see Appendix IB2.3K it is perhaps 
illuminating to provide an individual derivation for this quantity. 

3.5 Tension-Force Power Spectrum 

The tension force [Eq. 10] is F, = BjdjBj = djHjj, where we have 
omitted the factor of 1/4tt). Therefore, its spectrum [Eq. l[3}] is 

T{k) = 4xk 2 ®(k), (17) 
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where 



©(*) = <F*(A)F,(A)> = AXf.„.(A>. 



(18) 



We do not have any directly observable information about k z j= 0, 
so let us set A = k ± = fc ± (cos ip, sin tp, 0). Then 



fl>(A ± ) - <D] + <D 2 



(19) 



where Oi is the part that is directly recoverable from the Stokes 
correlators [using Eq. d 1 5 b 1 : 

®i = *2 |c..<A ) + <\. ...A .| + e y \(\ . (A ) + C I A )| 
+ 2k x kyRe |c,,,,(A> + C , (A.)| 

= - k 2 [E„ + Egg + E(/(7 + 2Re (E/g cos 2</> + E ;I / sin 2</>)] , 

(20) 

whereas 2 is the part that contains magnetic-field components 
parallel to the line of sight and, therefore, not picked up by the 
polarized-emission observations: 



$2 = KC w {kJ + ktCy^ikJ + 2k s k y ReC t: , yz (kJ. 



(21) 



y(*x) - C. (A )| 



It is to reconstruct this missing information that we have to assume 
isotropy, because it gives us a symmetry relationship between the 
unobservable correlators and the observable ones. If no mean field 
is present (B = 0), it is possible to show (see Appendix lB2,3b that, 
for a statistically isotropic magnetic-field distribution, 

*2 = ^{c w (fc x )--^-[c x 

= - k 2 (L UL i - T.QU tan 2tp) . (22) 

Assembling the directly observable [Eq. ( 120| >1 and the inferred 
[Eq. d22H part of the tension-force power spectrum, we arrive at an 
expression for O(fc) solely in terms of the Stokes correlators. There 
are two further steps that need to be taken to bring this expression 
into a practically computable form. 

Firstly, let us recall that, while the Stokes correlators in 
Eq. l |20| > and Eq. d22t depend on the vector k ± , the tension-force 
spectrum C> must depend only on k = \k±\. It is, therefore, permis- 
sible (and, in fact, increases the quality of the statistics) to average 
our expression for O over the angle tp (i.e., over a shell \k ± \ = k in 
the wavenumber space). 

Secondly, the fact that, for an isotropic field, only 4 of the 
6 available Stokes correlators are independent [see Eq. J16H can 
be used to construct many theoretically equivalent expressions for 
<b(k). Additional freedom comes from the angle independence of 
<t>(k) and, therefore, the possibility of doing weighted angle aver- 
ages (see Appendix IB2,2t . The strategy for choosing a particular 
formula for practical computations is to avoid having singularities 
in the coefficients: such as the factor of tan2</? in Eq. d22t . How to 
do this systematically is explained in Appendix IB2. 21 but here we 
simply give the result: 



T(k) 



5*f * 



£// + 2 (E/g cos 2tp + E/ u sin 2tp) - Eg;/ sin4</j 

(23) 



+ - (3 - cos 4<p) E ee + - (3 + cos Atp) X vu 



This formula is derived in Appendix lB 2 . 3 I from our general method, 
but can also be easily seen to follow directly from Eq. d20| l and 
Eq. d22t via Eq. dl6) , angle averaging and multiplication by the 
wave-number-space volume factor of 4nk 2 [see definition of T(k), 



Eq. d!7bl. Eq. d23t is our final expression for the tension-force 
power spectrum. 

Thus, we have accomplished our goal of showing that, despite 
the scarcity of the observable information, the tension-force power 
spectrum can be fully reconstructed from the available Stokes cor- 
relators (in Appendix |B2.2| we also show how to construct all other 
observable 4th-order quantities). In Sec. [4] we will test our method 
of doing this, but first, we generalize it slightly to the case of weak 
mean field. 



3.6 Case of Weak Mean Field: Observing the Magnetic-Field 
Power Spectrum 

We now relax the assumption that B = in Eq. (|5j- Then the 4th- 
order correlation tensor C,j_„„, [Eq. l l 1 I M can be written in terms of 
the mean field and of the correlation tensors of the fluctuating field: 

C,j, m „(r) = B^jB m B n + B-Bj{b' m b' n ) + B m B„(b i b j ) 
+ Btf m {bjb' n ) + Bfi n (bjb' m ) 
+ B J ^ m {b i b' n ) + BjB n (b j b' m ) 
+ Bi<b j b' m b' n ) + Bj{b i b' m b' n ) 
+ B m (bibjb' n ) + B n (bibjb' m ) 

+ {bibjb'M, (24) 

where unprimed quantities are evaluated at x and the primed ones at 
x + r. Due to homogeneity, correlation tensors depend only on r and 
not on x (and the statistical average can be interpreted as a volume 
average over x). This means that the first three terms in Eq. d24t 
have no spatial dependence at all, while the rest of the tensor can 
be written in Fourier space as follows: 



C ijjnii (A) 



= B^ m Cj, n (k) + B^ n Cj^k) 
+ B J : B m c ij ,{k) + BjB„c ij „(k) 
+ B i c* mntj (k) + B j cl n:i (k) 
+ B ni Cij n (A) + B„Cij, m (k) 

+ c '(7.imi(A). 

This is the Fourier-space correlation tensor introduced in Eq. J 1 2b . 
which has now been expressed in terms of the mean field and the 
second-, 3rd- and 4th-order correlation tensors of the fluctuating 
field: 



(25) 



c,,„,(A) 



,(*) 



(b*(k)bm(k)) 



(hUk)b m (k)) 



r (b,b'J, 



tfre-'^ibibjb'J, 



(26) 



(27) 



(28) 



c,;,,„„(A) = (h^kViUk)}, = ^ J ' d i re-' kr (b,b ] b' m b'„), 

where Ay (A) = (I IV) j ePx e- ikx b t {x)bj{x). 

Thus, the presence of the mean field leads to second- and 3rd- 
order statistics of the fluctuating field appearing alongside the 4th- 
order ones in the tensor Cjp„„. Since the Stokes correlators probe 
the total field, this means that some information about the second- 
and 3rd-order statistics could be extracted from them, provided the 
mean field itself can be independently determined and thus used 
as a "probe" (in fact, it turns out that only its orientation generally 
has to be known and even that knowledge is not always necessary, 
although easily obtainable; see Appendix |B3t . 

As before, we need additional symmetry assumptions about 
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the fluctuating field in order to make a transition from the Stokes 
correlators to theoretically/physically interesting quantities. The 
technically rigorous choice would be to assume that the statistics of 
b will depend on one special direction, that of the mean field, and 
be isotropic in the plane perpendicular to it. This, however, leads to 
a very large number of independent scalar functions appearing in 
the general form of C, ;> „„, and while it is probably worth working 
them all out, it is quite unlikely that the 6 available Stokes corre- 
lators will be sufficient to reconstruct anything of value. Therefore, 
we make a simplifying assumption (the physical grounds for which 
are discussed in Sec. 13.21 ) that the mean field is so weak (B <k (b 2 )) 
that the fluctuating field remains statistically isotropic. Under this 
assumption, the case of a weak mean field becomes a straightfor- 
ward generalization of the zero-mean-field case considered above. 
The main gain is that a weak mean field allows us to use the Stokes 
correlators to determine not just the power spectrum of the tension 
force but also the power spectrum of the magnetic field itself: as 
M(k) = Antfcij [Eq. {!}], it is recovered from the second-order 
terms in Eq. d25| l. 

The mathematical details of reconstructing the magnetic-field 
power spectrum are relegated to Appendix IB3.2I There are many 
equivalent expressions that can be derived for it; here we display 
three of them: 



M(k) 



2k 2 



—2 






cos 20 • 




2k 2 


—2 






sin 20 «• 




k 2 



f 

Jo 



dtp ReS;, 



2B sin 20 J o 



pin 

I dlfi (I.QQ - I.UU) ; 

Jo 



(29) 



where is the angle between the x axis and the projection of the 
mean field onto the plane perpendicular to the line of sight, B ± is 
the magnitude of this perpendicular projection. Although these are 
easy to measure (Appendix IB3. 3b . they are manifestly not neces- 
sary to determine the functional shape of the spectrum. Thus, we 
have three independent expressions from which we can deduce this 
functional shape. That the results should be consistent with one 
another is a good test of our assumptions (most importantly, the 
statistical isotropy of the fluctuating part of the field). 

The calculation of the tension-force power spectrum is entirely 
analogous to the zero-mean-field case (see Appendix IB3,5I >, The 
result is that Eq. j23\ still holds subject to two modifications: real 
part has to be taken of all Stokes correlators and a term proportional 
to M(k) has to be subtracted, namely 



T^ (k) = Re7^(Jfc) - - k 2 B 2 ± M(k) 



(30) 



where Tg =0 (k) is given by Eq. {23}. 

Finally, a disclaimer is in order with regard to the practical ap- 
plicability of the results obtained for the case of a weak mean field. 
Since we assumed the mean field to be small compared to the fluc- 
tuating field, B <K (b 2 ), the terms in Eq. i24l that contain 6, are 
small compared to {bibjb' m b' n ). Thus, in order for the second-order 
statistical information in Eq. d24l > to be recoverable, the errors as- 
sociated with the imperfect isotropy of the fluctuating field must be 
very small — smaller than 0(B ). It is not guaranteed that this is ei- 
ther justified physically or achievable in practice and the verdict on 
the usefulness of the results of this section will depend on extensive 
numerical tests, which will not be undertaken in this paper and are 
left for future work. 



4 NUMERICAL TESTS 

Having presented the analytical derivation of our method, we now 
present a proof-of-concept numerical test by analyzing two data 
cubes containing randomly tangled magnetic fields: a saturated 
magnetic field ge nerated by fluctuation dyn amo in an MHD simu- 
lation (Run S4 of lSchekochihin et alj|2004> and a divergence-free, 
random-phased Gaussian field synthetically generated to have the 
same spectrum as the MHD field (Fig.|2j snapshots of the two fields 
are shown in Fig.QJ. Both fields have zero mean, so the results of 
Sec.|3.6|are not tested here. 



4.1 Case of p = 3 

We first test the validity of our method for the case of the electron 
spectral index p = 3, assumed throughout the analytical develop- 
ments presented above. For each data cube, we designate one of 
its axes as the "line of sight" and integrate the field along it ac- 
cording to Eq. ((9}. This produces three synthetic two-dimensional 
Stokes maps (Fig.[3j examples of such /, Q and U maps are shown 
in Fig. [4}. Since we have the full three-dimensional information 
for both fields, we can compute the tension-force power spectra di- 
rectly and then compare them to the spectra obtained by applying 
our estimator, Eq. {23}, to the synthetic Stokes maps. 

In Fig. [5] (left panel) we plot the tension-force power spectra 
reconstructed using our estimator, Eq. {23}, for a realization of an 
MHD simulated field and for a synthetic Gaussian field. They are 
compared to the same spectra directly computed from the three- 
dimensional data cubes [according to Eq. {3}]- The reconstructed 
spectra are obtained as an average over three synthetic Stokes maps, 
each obtained by choosing as the "line of sight" one of the three 
orthogonal axes of the data cube. This allows us to estimate the 
accuracy of the reconstruction, represented in Fig. [5] by the error 
bars. 

For both types of field, the performance of our estimator is 
clearly excellent. The relative error bars for the Gaussian random 
field are substantially smaller than for the MHD field, which makes 
sense in view of the former's more small-scale and less structured 
character. The salient point that emerges from the comparison of 
the two fields is that the tension-force spectrum can be recovered 
from the synthetic observations with an accuracy easily permitting 
to discriminate between the structured (folded) MHD field and the 
structureless Gaussian one. This suggests that the proposed estima- 
tor is a robust tool for diagnosing magnetic turbulence from po- 
larized emission data and for discriminating between different sce- 
narios of magnetic-field evolution and saturation (see discussion in 
Sec.|2}. 

The test that we have presented only allows us to assess the 
quality of our method under idealized conditions, namely, assum- 
ing that the observation is noiseless, that no observational-window 
effects are present, that the relativistic-electron energy distribution 
is homogeneous and has the spectral index p = 3, and that the Fara- 
day rotation is either negligible or has been effectively subtracted 
(Sec. l3.3l t. Thus, the errors in our reconstructed spectra are due to 
two factors: firstly, a certain amount of information is lost in the 
projection of the three-dimensional data onto a two-dimensional 
Stokes map (the line-of-sight integration); secondly, the assump- 
tions of statistical homogeneity and isotropy (Sec. |3.2K upon which 
our estimator depends, are imperfectly satisfied by any particular 
realization of the field. It is interesting to ask what is the relative 
contribution of these two sources of inaccuracy to the errors of re- 
construction represented by the error bars in the left panel of Fig. [5] 
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MHD Synthetic Gaussian 




U, p = 3 U, p = 3 

Figure 4. Examples of synthetic /, Q and U maps generated from an MHD-simulated magnetic field (left panels) and a synthetic Gaussian field (right panels). 
The same data cubes were used as in Fig.[T] 
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k/2TT k/2TT 

Figure 5. Left panel: The bold solid red line shows the tensio n-force power spectrum rec onstructed via Eq. (23) from the synthetic Stokes maps (Fig. [4) based 
on an MHD-simulated field (saturated state of the Run S4 of Schekochihin et al. 2004]). The bold dotted black line is the same spectrum computed directly 
from the full three-dimensional data (same as Fig. [2] right panel). The errors bars on the estimated spectrum are obtained by comparing results from synthetic 
Stokes maps obtained by integrating along three orthogonal "lines of sight" (the three axes of the data cube). The thin solid blue line with error bars and the 
thin dotted black line represent analogous information for a synthetic Gaussian field. Right panel: Similar to the left panel, but the reconstructed tension-force 
spectra are based not on the estimate Eq. 423t but on the full information about the projected (line-of-sight integrated) spectra, i.e., they are given by the sum 
of® i [Eq. I20H and ®2 [Eq. 121 H calculated in terms of C,-y, mn (^±) (including its unobservable line-of-sight components). 



This is addressed the right panel of the same figure, which is anal- 
ogous to the left panel, but instead of the spectra reconstructed via 
Eq. ( 123b . it shows the spectra resulting just from the line-of-sight 
integration (setting k z = 0) of the full tensor C, ;.,„„, i.e, they use 
the unobservable z components of this tensor that enter in Eq. d21t 
rather than infer them from the observable components and the 
isotropy assumption. Comparing the right and left panels of Fig. [5] 
suggests that much of the reconstruction error (especially at large 
wave numbers) is due to the loss of information associated with the 
line-of-sight integration, not to imperfect isotropy — and this is de- 
spite the fact that the MHD field contains magnetic structures with 
virtually box-size parallel coherence lengths (see the left panel of 
Fig.[Q- 



4.2 Case of p t 3 

Assuming that the electron spectral index p = 3 was an idealization 
of the real observational situation that we needed for the theoretical 
justification of our method because Stokes parameters are strictly 
quadratic in the magnetic field only if p = 3 (see Appendix lAl >. 
While p = 3 is not a bad approximation of reality, one cannot ex- 
pect it to be satisfied very precisely (see discussion and references 
in Sec. 13. 3I >. so in order for our method to be practically useful 
for real observations, it must be reasonably insensitive to the exact 
value of p. This sensitivity is very easy to test. 

Let us generalize our definition of the Stokes parameters 
[Eq. l©] to the case of p + 3: suppressing the dimensional pref- 
actors as before, we get (see Appendix [At 

/(xj = 1 J dz[Bl(x) + Bl(x)f'~ m [Bl(x) + Bl(x)], 
2(Xx) = i J dz[Bl(x) + B 2 y (x)f'~ m [Bl(x)-Bj(x)], 

U(x ± ) = dz[Bl(x) + B 2 y (x)f~ 3m 2B x (x)B y (x). (31) 

Clearly, for p > 3, the extra factor of (B\ + flJ)(p-3)/4 causes me 
statistics to be effectively weighted towards regions where the field 



is stronger, for p < 3, towards those where it is weaker. This point 
is illustrated by Fig. [6] which shows that increasing/decreasing p 
roughly corresponds to increasing/decreasing the contrast in the 
Stokes maps. 

The range of values that p can realistically be expected to 
take is roughly p e [1.5,3.5] (see references in Sec. 13.3b . Since 
this implies that (p - 3)/4 e [-0.375,0.125] are not very large 
powers, there is a priori a hope that the effect of deviations from 
p = 3 might not be catastrophic for our estimator. This, indeed, 
proves to be correct. In Fig. [7] we show the tension-force spectra 
reconstructed from Stokes maps generated using Eq. d31b with a 
number of values of p and compare them to the true spectra. Even 
for values of p significantly different from 3 (roughly in the range 
p e [2.5, 3.5]), our estimator works extremely well, except at the 
highest wave numbers. 

Note that the extra factor of (B 2 X + 5j) (p_3)/4 in Eq. ([3B changes 
the overall amplitude of the Stokes parameters in comparison to 
what it would have been with p = 3, so we can only hope to 
recover the functional shape of the tension-force power spectrum, 
not its overall magnitude. In the numerical data used above this 
potential source of reconstruction error is not very visible because 
values of the magnetic field are close to unity in code units, but 
in any realistic observational situation, the shift in amplitude of 
the Stokes parameters may be significant. Importantly, however, 
we see in Fig. [7J that in all cases we have tested, the shape of 
the reconstructed tension-force power spectrum still makes it 
unambiguously possible to discriminate between qualitatively 
different field structures as represented by the MHD and Gaussian 
fields. 

The numerical tests presented above are meant to demonstrate 
in principle that the approach taken in this paper is a valid one. We 
did not attempt to test the robustness of our approach by including 
into our synthetic data model all of the complications that will arise 
in handling real observational data. A known caveat is that obser- 
vational window functions, due to the finite size of the radio source 
of the telescope beam, will lead to a redistribution of power in the 
recovered spectrum, so that the large-scale power may swamp the 
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U,p = 1.5 U,p = 4.5 

Figure 6. The Stokes maps calculated according to Eq. j3i\ for two extreme values of the electron spectral index, p = 1.5 (left panels) and p = 4.5 (right 
panels). These are to be compared with the Stokes maps for p = 3 shown in the left panel of Fig. [4] 
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Figure 7. These plots show the same comparisons as the left panel of Fig. [5] but for a number of values of the electron spectral index p ?= 3. The Stokes maps 
were calculated according to Eq. Some of these Stokes maps are shown in Fig. [6] 
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signal at high wave numbers dVogt & EnBlinll2003h . Both further 
tests and applications of our method to real data are left to future 
work. 



5 CONCLUSION 

We have demonstrated that it is possible to reconstruct the power 
spectrum of the tension force associated with tangled astrophysical 
magnetic fields as a linear combination of the radio synchrotron 
observables, the Stokes correlators. This was done under a set 
of simplifying assumptions about the synchrotron emission data 
(Sec. |3.3b and also by assuming a statistically homogeneous and 
isotropic stochastic magnetic field (Sec. 13.21 ). The tension-force 
power spectrum emerges as a particular case from a subset of ob- 
servable 4th-order statistics (Sec. |3.5| and Appendix |B2.3| l — a non- 
trivial fact because in general, the Stokes maps do not carry suffi- 
cient information to reconstruct all of the 4th-order correlators of 
the magnetic field (see Sec. |3.4| and Appendix |B2.2t . 

The observability of the tension-force power spectrum is a 
stroke of good fortune because this quantity plays an important 
role in diagnosing the spa tial structure of the magnetic turbulence 
dSchekochihin et al]|2004h and allows one to distinguish between 
different theoretical scenarios for the evolution and saturation of 
the cosmic magnetic field, which was not possible to do on the ba- 
sis of lower-order statistics such as the magnetic power spectrum; 
it also reveals physically interpretable dynamical properties of the 
system under observation, namely the force exerted by the field on 
the ambient plasma (see discussion in Sec.|2j. 

Furthermore, we have shown that if the observed magnetic 
field possesses a small regular component that does not affect the 
isotropy of the fluctuating part of the field, it may be possible to 
obtain from the Stokes maps the power spectrum of the fluctuating 
field itself, as well as that of its tension force (Sec. 13. 61 . 

Thus, physically relevant information about the spatial struc- 
ture and dynamical properties of the magnetic turbulence is con- 
tained in the polarized emission maps and can be extracted. This 
work is an attempt to pave the way towards analyzing the large 
amount of exi sting and upcoming radio- synchrotron observa tional 
data (see, e.g., Gaensler 2006l lEnfllin et al.ll200dlBec3l2008l) with 
the aim of achieving a better understanding of the nature of magne- 
tized turbulence in cosmic plasmas. 
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APPENDIX A: SYNCHROTRON EMISSION AND THE 
STOKES PARAMETERS 

A spatially homogeneous, pitch-angle-isotropic and power-law 
distributed in energy relativistic-electron population is assumed 
[Eq. ©], The resulting synchrotron emission is partially linearly 
polarized. Its intensity and polarization depend solely on the mag- 
nitude and orientation of the magnetic field B ± projected onto the 
plane perpendicular to the line of sight and on the electron distri- 
bution [Eq. ©]. 

The synchrotron emissivity (i.e., power per unit volume per 
frequency per solid angle) is usually subdivided into two com- 
ponents, respectively perpe ndicular and parallel to B ± : following 
iRvbicki & LightmaiJ ( ll979h . 



j ± (co,x) 



[F(p) + G(p)}w {l -< m \BAxi p+im 
[F(p)-G(p)] co fi -"> l2 \B ± (x)\ { " +ly2 , 



(Al) 



where oj = 2nv, v is the observation frequency, x is the spatial 
position, p is the spectral index of the electron distribution [Eq. ([8}] 
and 



F(p) 



G(p) 



V3e 3 ( 2m P c 



d-p)/2 



C 



32n 2 m e c 2 \ 3e 



V3* 



12/ p+ 1 



2m P c 



(1-/0/2 



c 



32n 1 m e c 1 \ 3e 

( P Wi'-^ 2 T\'- + 



x r 



12 



12 



(A2) 



where m e is the electron mass, e is its charge, c is the speed of light, 
and C is the prefactor of the electron distribution [Eq. {§}]. 

The specific intensity I and the polarized specific intensity PI 



are given by the following line-of-sight integrals (see lBurnll 19661) : 

J "here 
dz[j\(co,x) + j\(ai,x)], 
there 



J there 
-here 



PI(w. 



/-»nere 

x ± ) = dz[M<o,x)-j ll (a>,x)]e- 2i x<*\. 

<J there 



(A3) 



where z is the line-of-sight coordinate, x x = (-*, y) is the position 
vector in the plane of the sky (perpendicular to the line of sight) 
and the polarization angle is given by 



X (x) = xo(x) + A 2 RM(x\ 



(A4) 



where A = c/v is wavelength of the observed emission, the intrinsic 
polarization angle is 



_1 By 

Xo(x) = tan — , 

and the Faraday rotation measure 



RM(x) 



2mn 2 ,c 4 J tha 



Ehei 



dzn the (x)B-(x), 



(A5) 



(A6) 



where n^ e is the density of thermal electrons and B : the projection 
of the magnetic field on the line of sight. 

The Stokes parameters are now defined as follows 



/ = dm, Q-iU 



I 



rfDPI, 



(A7) 



where the integration is over the solid angle of the angular resolu- 
tion element (the observational beam). If the spectral index is taken 
to be p = 3 and the Faraday rotation in Eq. dA4l is assumed to be 
negligible (as discussed in Sec. l3.3t . the Stokes parameters depend 
quadratically on the components of the magnetic field perpendicu- 
lar to the line of sight, B x = B ± cos^ and B y = B ± sin^- - Indeed, 
using the above definitions, we get 



//->here 
dQ. I dz{B 2 x + B 2 ), 
*J there 

//^here 
dQ. I dz(B 2 x -B 2 ), 
■J there 

fdnf 

J Jlhe 



2G(3) uf 



dz2B x B y . 



(A8) 



From these formulae, we recover the analytically convenient defi- 
nitions of the Stokes parameters, Eq. j5J, by dropping the dimen- 
sional prefactors and the integration over the angular resolution el- 
ement and normalizing the integrals by the depth of the emission 
region. 

Somewhat more generally, for arbitrary p, but still neglecting 
the Faraday rotation, we have 

/ = 2F(p) oj"-"" 2 f dQ f hCre d Z (Bl + B^f" m (Bl + B 2 ) , 

*J <J there 

Q = 2G(p) W (1 -"" 2 [dQ r e dz{Bl + B 2 ) ( "- 3)l4 (Bl-B 2 ), 

-J *J there 

U = 2G(p)co (, -" )l2 f dQ f m dz(B 2 x + B 2 ) { "^ VA 2B x B y . 

*J <J there 



(A9) 



These formulae are the basis for Eq. J3 lb . 
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APPENDIX B: FOURTH-ORDER CORRELATION 
TENSOR AND ITS REPRESENTATION IN TERMS OF 
STOKES CORRELATORS 

In this Appendix, we derive the general form of the 4th-order cor- 
relation tensor C,,.,,,,, [Eq. i ll 1H for a statistically homogeneous and 
isotropic magnetic field and show what part of the relevant statisti- 
cal information can be recovered using Stokes correlators. 



Bl Symmetries and the General Form of C, ;m „ 

In Eq. l !25t . the tensor C,j„,„ is written in Fourier space in terms of 
the mean field B, and of the second-, 3rd-, and 4th-order correlation 
tensors of the fluctuating field b h denoted c j m , Cy iW and Cy iMB . Each 
of these correlation tensors de pends on a certain number of scalar 
correlation functions (see, e.g.. lRobertson|[l940l) . This number can 
be constrained if we take into account some intrinsic properties of 
correlation tensors (permutation of indices), of the field they are 
constructed from (it is a real, divergence-free field), and additional 
symmetries we assume (homogeneity and isotropy). Let us imple- 
ment these constraints. The procedure is least cumbersome when 
applied to the second-order correlation tensor. We will explain it in 
detail on this example and then proceed analogously with the 3rd- 
and 4th-order correlators. All further calculations will be in Fourier 
space, but exactly analogous calculations can be done in position 
space if it is necessary to compute position-space correlators. 



Bl.l Second-Order Correlation Tensor 

For a statistically isotropic field, the second-order correlation ten- 
sor depends on three scalar functions — this is shown by construct- 
ing c,,,„ out of all possible isotropic second-rank tensors. In three 
dimensions, the available building blocks for these tensors are (5„„, 
e imp and k h the unit vector in the direction of k. Therefore, 

c.>(*0 = j \ m i 6 im + m 2 kikm\ + im 3 £ imp kp, (Bl) 

where the scalar coefficients m\, m 2 , m 3 can only depend on k = \k\. 

Since bj(k) is a Fourier transform of a real function, we must 
have bj(-k) = b*(k), whence 

cU-k) = c* m (fc). (B2) 

It is easy to see that this implies that m\, mj and m 3 are real (the 
factor of i in front of 013 was chosen deliberately to arrange for this 
outcome). 

Since c, „, is a correlation tensor, it has a symmetry with re- 
spect to permutation of its indices: 

<W(*) = (b'Jk)bi(k)) = <„,(*). (B3) 

This does not bring any new information beyond the reality of mi, 
m 2 and m$ . 

Finally, the magnetic field is solenoidal, kjbj(k) = 0, so we 
must have 

kjCj m = k m Cj m = 0. (B4) 

This gives mi = —mi, so the general form of the second-order cor- 
relation tensor is 

Ci,,„(k) = - mi(k) (dm - kik,„) + im 3 (k)€ imp k p , (B5) 
i.e., it depends only on two scalar functions. If we take the trace 
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of this tensor, we obtain the magnetic-energy power spectrum 
[Eq. QJ]: 

M(k) = 4mt 2 c;.,(/t) = 4nk 2 mi (k), (B6) 

so we do not need to know m 3 if we are only interested in the power 
spectrum. IVogt & Enfilml J20031. I200I) used this property to pro- 
pose a way to measure the magnetic power spectrum solely in terms 
of the scalar correlation function of the Faraday rotation measure 
associated with a given magnetic-field distribution: although only 
one scalar function was available this way, assuming isotropy and 
restricting one's attention to a particular quantity of physical inter- 
est made it possible to make do with incomplete information. We 
follow the same basic philosophy in this paper, primarily as applied 
to the 4th-order statistics. 

Note that is a measure of reflection (parity, or mirror) non- 
invariance of the magnetic field. If OT3 ^ 0, the field has helicity. If 
we demand mirror symmetry of the field, 

c,, m (-k) = c Um {k), (B7) 

we find 7«3 = 0. We will see that normally we do not have to make 
this assumption because in many cases, the mirror-noninvariant 
terms are not present in the quantities of interest (as was the case 
with the power spectrum). 

B1.2 Third-Order Correlation Tensor 

Analogously to the above, we construct the general isotropic 3rd- 
order tensor as follows 

(k) = i( 

-\-{Z(~,€jjpkpk m + ci-j€j m pkpkj + a$€j m pkpkj, (B8) 

where a, , . . . , a 8 are functions of k = \k\ only. 
Reality of the fields h is and b„, implies 

cnJrK = c*j, m (k), (B9) 

whence ai, . . . , as are all real. 
Permutation symmetry, 

Cji, m (k) = c m (k), (BIO) 

implies 02 = a 3 , fls = fl6 = 0, and a-] = a$. 

Solenoidality of the magnetic field implies 

k m c iJtm (k) = 0, (Bll) 

whence ai = and a 4 = -2a 2 . 

Thus, the general form of the 3rd-order correlation tensor is 

Cijjnik) = ia 2 (k)(s im kj +6 Jm k' i -2k' i k'jk^ 

+ a 7 (k) (e imp k p kj + ej mp % p k\j . (B12) 

B1.3 Fourth-Order Correlation Tensor 

In the 4th order, the number of terms in the general tensor becomes 
quite large. Constructing this general form out of the usual building 
blocks, Sij, k, and e !jm , we get, via straightforward combinatorics, 

(k) = 

+ c^8ijk m k n + c$S nm kjkj 

"I" Csdimkjkn ^'l^iiikjk,,, + Cgdjmkjk,, + CgSjnkjk,), 

+ C\Qkikjk m k n 

+ i ( C[i €jj m k )} + C\2^ijnk m + C[2^imnkj + C\4€j mn kj 
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+ c lS^inpkpSj m + C 196 JmpkpS j n + C2Q€jnpkpSj„, 

+ C 21 €ijpkpk m kfj + C22^mnpkpkjkj + C23^impkpkjk„ 

+ ^24^inpkpkjk m + C25^jmpkpkjk n + C26 Ejnp^p^k i^m J j 

(B13) 

where c\, . . . , c 2 6 are functions of k = |fc| only. Note that there are 



no terms of the form £ij P k p e„ mq k q because 

£ijp£mnq = S im 6 )„6pq + 6j„S jq6 pm + 6 iq 6 jmd ,,„ 
- 8j m Sj q 6 pn - 6 in 6j m 6pq - S iq Sj„6 pm , 



(B14) 



so such terms are already present in the general form we have con- 
structed. 

Reality of the field hjj = b,bj implies 



Cij,mn{ k) — Cjj mn {k\ 

whence ci, . . . , C26 are all real. 

There are three permutation symmetries: 



gives c 2 = c 3 , c 6 = c 8 , c-i = c 9 , c u = Cu = c l5 = c 2 i = 0, c l3 

c n = C[9, c 18 = C-2Q, t'23 = c 25, C 2 4 = C26, 



(B15) 
(B16) 



L 1 /.inn 



(*). 



gives additionally c 6 = c 7 , c g = eg, c 13 = c 14 = c 16 
C17 = Ci 8 , C19 = c 2 o, t'23 = t'24, c 2 5 = c 2 6, and, finally, 



(B17) 

c'22 = 0, 

(B18) 



gives c 4 = c 5 . 

Assembling all this information, we find that the general 4th- 
order correlation tensor only depends on 7 scalar functions: 

c ijjrm(k) = c l(k)SjjS„,„ + C 2 

+ c 4 (k) (Sjjk,„k„ + S„ m kik^ 

i m kjk n + 5j n kjk m + 5j in kik fl + Sj n kjk m ^ 
+ c [0 (k)ki kjkmkn 
+ icn{k)\etmpkpdj n + e inp k p Sj m 

£ji>i P k p Sjn + £jn P k p Sjt]i\ 



+ ic^iiji) i^im p k p kjk n + €{ np k p kjk m 



~t~ ^jmpkpkik n + £ji, p k p kjk m ^j . 

B2 Observables in the Case of Zero Mean Field 



(B19) 



Let us first examine the case B = 0, so we are only concerned 
with the 4th-order statistics. We will need explicit expressions for 
the coordinate-dependent components of the tensor c, Jim „ in terms 
of the coordinate-invariant functions ci, c 2 , c 4 , C6, cio, cn and C23 
[Eq. ( IB 1 9I >1 , As the polarized emission data on the magnetic field 
arrives in the form of line-of-sight integrals (Sec. I3.4t . we have to 
set k z = everywhere — no information on the field variation in this 
direction is available. However, because of the assumed isotropy, 
the dependence of the invariant scalar functions on k ± = \k ± \ con- 
tains the same information as their dependence on k = \k\. Let us 
denote by <p the angle between k ± and the x axis. This means that 
we set k = (cos tp, sin tp, 0). Then the components perpendicular to 
the line of sight are 

c xx ,xx(k ± ) = ci + 2c 2 + (2c 4 + 4c 6 ) cos 2 tp + c w cos 4 tp, 



c yy,yy(k±) 

c xy,xy(k±) 
c xx^\y(k±) 
c yy r xy(k±) 



Ci + 2c 2 + (2c 4 + 4c'6) sin" <p + cio sin <p, 

C\ + c 4 + Cio sin 2 <^> cos 2 

. 7 2 
c 2 + eg + c^o sin" tp cos" 

(c 4 + 2c 6 ) cos ^ sin^) + cio cos 3 <p sin<p, 

(c 4 + 2c 6 ) cos </) sin </) + cio cos ip sin 3 1/3. 



(B20) 



These are the only components of c, ; „,„ that are directly sampled by 
the polarized emission. The components parallel to the line of sight 



arc 




r- 


= 


Ci + 2c 2 , 


^ XX 


= 


C[ + c 4 cos 2 ^1, 


Cyy 


= 


ci + c 4 sin 2 


Cxz 


= 


2 

C2 + C6 COS <p, 


Cyz 




c 2 + c'6 sin 2 


Cxz 


= 


eg sin 1/3 cos tp. 



(B21) 

Information about these components can only be obtained by rely- 
ing on the isotropy assumption as they are expressed in terms of the 
same invariant scalar functions as the perpendicular components. 
Note that setting k- = has led to all information being lost about 
the mirror-asymmetric part of the tensor, so no quantity involving 
Cn or C23 can ever be reconstructed from polarized emission. 



B2.1 Stokes Correlators 



Using Eq. JB201 > and the expressions for the Stokes correlators 
given by Eq. l !14b . we get 



Zuu(k ± ) 
S/ fl (*j.) 



4 (Ci + Ci + C 4 + C6) + Cio, 

4 (c 2 + C(,) + cio cos 2 2^, 
4(c 2 + c 6 ) + C|osin 2 2y3, 
(2c 4 + 4c6 + cio) cos 2tp, 
(2c 4 + 4c6 + cio) sin 2^, 
Cm sin 2(/>cos 2y, 



(B22) 



Note that Hjq and Z/y contain the same information and so do Hqu 
and Hqq - The relations between them follow immediately 
from Eq. dB22b and are given by Eq. 1 II6I 1. 

Thus, only 4 of the Stokes correlators are independent: S/;, 
two of ~Lqq, Ey;/, and one of £/c/. We see that inEq. ( IB 22t . 
these 4 independent observables are expressed in terms of 5 invari- 
ant scalar functions, ci, C2, c 4 , C(, and cio, which cannot, therefore, 
all be reconstructed from polarized emission data even if isotropy is 
assumed (as we explained earlier, two other functions, c 17 and c 2 3, 
which complete the full set of 7 alluded to in Sec. 13. 41 can never be 
known from polarized-emission data). 

The relations between the Stokes correlators and the invariant 
scalar functions given by Eq. ( IB22t contain angular dependence. It 
will be convenient for practical calculations to express all observ- 
ables in terms of angle averages (i.e., averages over all orientations 
of k ± ): 

Zi (it) 

E 2 (fc) 

£,(*) 

£*(*) 



■4(d + c 2 + c 4 + c 6 ) + cm, 

> (C'2 + C 6 ) + Cio, 



— J #(S Qe +S TO ) = 8( 

^- J"rfi^[(E ee - E t/ (/)cos4v3 + 2S e c/ sin4(^] = c 10 , 

— f dtp(l.]Q cos 2</j + Z/y sin 2^3) = 2c 4 + 4c'6 + c^o. 
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(B23) 

These formulae again give us 4 independent observable scalar func- 
tions, but now the quality of the statistics should be improved by the 
angle averaging. We will find that it is most convenient for practical 
calculations to use 2^, . . . , £4 as the basic set of observables (see 

Sec. nml 



B2.2 General Observable Quantities 

Thus, only some 4th-order statistical quantities are observable. It is 
not hard to work out the condition for them to be so. First, as we 
already explained in our discussion of Eq. dB20b . Stokes correlators 
carry no information about anything that involves the measures of 
mirror asymmetry of the field c tl and C13. Let us then restrict our 
attention to quantities that contain only the remaining 5 invariant 
scalar functions that determine the 4th-order two-point statistics 
[Eq. dB19b l. In general, we would be looking for a scalar function 
that has the form 



*(*) = fid(k) + f 2 c 2 (k) + f 4 c 4 (k) + f 6 c 6 (k) + /iodo(fc), 



(B24) 



where fi, ...,fio are some known coefficients, which can be func- 
tions of k. Let us try to express this quantity in terms of Stokes 
correlators: this amounts to finding coefficients a, fi, 7, <5, which 
can be functions of k and ip, such that 



<S>(k) = \ [o&n + P^qq + J^uu + S 

4 \ cos 2tp 



(B25) 



Using Eq. dB22b . we get 

$>(k) = ac\ + (a + B + y)c 2 + (a + ^jc 4 + (a + B + y + S)c 6 

+ — (a + fi cos 2 2p + y sin 2 2<p + 5) c [0 . (B26) 
Comparing this with Eq. dB24b . we get 



a certain redundancy of information available from the Stokes cor- 
relators, as expressed by Eq. J 1 6b . mean that there are, in general, 
many different ways of reconstructing observable quantities. In the- 
ory they are all equivalent, but in practice one has to choose one 
with the aim of reducing noise and offsetting the potentially detri- 
mental effect of singularities in the coefficients associated with fac- 
tors of 1 / cos 2ip and 1 / cos 4<p. 

One method, which we have found to be quite effective, of 
avoiding this problem is to pick as our basic set of 4 observable 
scalar functions not the Stokes correlators themselves but the com- 
binations of their angle averages given by Eq. dB23b , Repeating the 
procedure we have just followed, we seek 0(£) in the form 



$>(k) = ^ + o-iXi + 0-323 + 0:4X4) , 



(B30) 



where the coefficients a { , . . . , a 4 are now functions of k only [there 
is no longer any angular dependence on either side of Eq. JB30H . 
Using Eq. dB23b . this becomes 

<i>{k) = ct\C\ + {a\ + 2a 2 )c2 + |oi + "^~j c 4 + ( a i + 2c*2 + o"4)c6 

+ -((*] +o 2 +03 +o 4 )c 10 . (B31) 
4 

Comparing this expression with Eq. dB24t as before, we get 



«i = fu 

ai = 5(6-/1). 

03 = j (3/i -h -4/4 + 8/ 10 ), 

0-4 = 2(/ 4 -/i) = / 6 -/ 2 , 



(B32) 



and the observability criterion is again given by Eq. JB28b . 
Eq. JB30t together with Eq. dB23b and Eq. dB32t give another ex- 
pression for a general observable 0(£), defined by Eq. ( IB24t and 
subject to the constraint Eq. l lB28t . 



fu 

l (f f 3/i -A -4/4 + 8/10 
2 I/2-/1 + - 



y 

6 



2(/ 4 - 



■ fi- 
ll) 



cos4<p 

3/i -A- 4/4 + 8/10 



cos4i/> 



fe-fi- 



(B27) 



The last formula gives two expressions for 6. In general, they do not 
have to be compatible and if they are not, the quantity ©(£) cannot 
be expressed in terms of Stokes correlators. Thus, we have derived 
a simple criterion: only those quantities given by Eq. dB24t 
are observable for which 



0. 



(B28) 



2/i-/ 2 -2/ 4 +/ 6 

If Eq. lB28t is satisfied, Eq. dB27t and Eq. {B25) give us 
a practical method for calculating O. As any interesting physical 
quantity O has to be independent of the angle ip between the wave 
vector k ± and the x axis of the coordinate system in which the 
Stokes parameters are measured, we are allowed to average over tp: 



®(k) 



1 

8^ 



dpW(p) 



aLu +fi^QQ + y^uu + S 



■-iQ 



cos Txp 



(B29) 



where W(ip) is some weight function, which must satisfy 
(1/2tt) £" dpW(p) = 1. The weighting does not theoretically af- 
fect the result, so W(p) can be chosen arbitrarily. 

The possibility of angle averaging with a weight function and 



B2.3 Tension-Force Power Spectrum 



In Sec. 13.51 we split the tension-force power spectrum into the 
directly observable part [Eq. l|20H, which could be recovered 
from the Stokes correlators without any assumptions, and the non- 
directly-observable part ©2 [Eq. i2U], which could only be recon- 
structed using some assumed symmetries of the 4th-order correla- 
tion tensor. Assuming isotropy, we infer from Eq. JB21t 



$2 = Ir(c2 + c 6 ). 



(B33) 



Using Eq. l |B20t to express c 2 + c b in terms of the perpendicular 
components of the tensor c,^ m „, we immediately recover Eq. d22t 
and the rest follows as explained in Sec. 13.51 

This was an ad hoc derivation specific to the tension-force 
power spectrum. Let us now demonstrate how the general method 
laid out in Sec. IB2.21 works for this quantity. 

Substituting Eq. JB 1 91 > into Eq. i fTSl , we get 



*(*) = k 1 [cdk) + 4c 2 (k) + 2c 4 (k) + 6c 6 (k) + c l0 (,kj] , 



(B34) 



a particular case of Eq. dB24t - The observability criterion given by 
Eq. dB28b is satisfied, so, using Eq. dB25t and Eq. dB27b . we obtain 



1* 



1 / 1 

Z//+TT 3 — 

2 \ cos 4ip 



2 \ cos 4ip 



2 

cos 2p 



-IQ 



(B35) 
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which can be angle-averaged with some weight function according 
to Eq. {B29} . 

An alternative expression in terms of averaged Stokes correla- 
tors follows from Eq. dB30l > and Eq. dB32| >: 



9(k) = \k 2 fa + ^ 2 - ±Z 3 + 2Z 4 



(B36) 



Substituting for £[ , . . . , E 4 from Eq. dB23b , we arrive at 

A' 2 



8?r Jo 
1 



Xji +2(Ljq cos 2tp + T. w sin 2tp) - E e( y sin 4tp 
1 



+ - (3 - cos 4tp) T.qq + - (3 + cos 4tp) ~L VV 



(B37) 



Our final formula for the tension-force power spectrum, Eq. {23}, 
follows from Eq. l |B37| > upon multiplication by 4nk 2 (the wave- 
vector-space volume factor). Note that the integrand in Eq. dB37t 
reduces back to Eq. dB35t if we make use of Eq. d 1 6b . but the ad- 
vantage of Eq. dB37| l is that it does not contain any singular coeffi- 
cients. 



B3 Observables in the Case of Weak Mean Field 

If a weak mean field is present, we proceed analogously to Sec. lB2l 
It is understood that the mean field is sufficiently weak so as not to 
break the isotropy of the fluctuating part of the field. Then, using 
Eq. {25}, Eq. dBU and Eq. dB12l > and setting k, = to express the 
line-of-sight integrals, we find 

—2 2 

C xx , xx {k ± ) = c xvx (k ± ) + 2B x sm tpm u 
< . I* I = c„,„.(k ± ) + 2b\, cos 2 ipm 1 , 
C xx , yy (k ± ) = c xx _ yy (k ± ) - 2B x B y sin tp cos tp m { 

+ 4i iyB x sin tp + B y cos tpj sin <p cos ip a^, 

1 I- - \ 2 

(\ (A > = c xytXy (k ± )+ -\B x cos<p-B y sincp) m\, 

C XXiXy (k ± ) = c XXiXy — B x sin (p \B X cos tp — B y sin tp) m i 
- 2i [b x sin tp - B y cos tp) sin 2 tp ai, 

Cyy, X y (k ± ) = c; vv _ , v (k±) + By cos tp \B X cos tp - By sin tp) m i 

+ 2i (b x sin tp - B y cos ip^j cos 2 tp 02, (B38) 



where the components of c,j >m „(A ± ) are given by Eq. dB20t and 
Eq. dB2"T} . 



B3.1 Stokes Correlators 



Using Eq. {B38} and Eq. d!4t . we find that the Stokes correla- 
tors are 



= 2 (b x sin tp - By cos tp^j my + 4th order, 
Y,QQ(k ± ) = 2 (b x sin tp + B y cos tpj mi + 4th order, 
J^uu(k±) = 2 (b x cos tp - By sin tp\ my + 4th order, 
~EiQ(k±) = 2 |s v sin 2 tp - B y cos 2 tp^mi 

- 4i [b x sin tp — B y cos tp^j sin 2tp 02 + 4th order, 
S /c/ (jtj_) = -|^ + B^sin2 V -2B v S v |m 1 

+ 4i [b x sin tp — B y cos ^ cos 2tp 02 + 4th order, 



^ou(k ± ) = - ^b] - B^j sin 2tp + 2B x B y cos 2tp^ m x 

- 4i(B x sin tp - By cos tp^j a 2 + 4th order, (B39) 

where the 4th-order parts of the correlators are given by Eq. {B22}. 

Thus, the Stokes correlators now contain not just the 4th-order 
statistics but also some information about the second- and 3rd-order 
correlation functions of the field, namely, the magnetic-field power 
spectrum m\ (k) [see Eq. dBlH and the 3rd-order correlation func- 
tion a 2 (k) [Eq. JB12H . We are not particularly interested in a 2 and 
notice that all terms containing it can be eliminated from Eq. dB39t 
simply by taking the real part of the Stokes correlators E/g, 
and ~Lqu- We will now isolate the second- and 4th-order contribu- 
tions to the Stokes correlators and calculate the power spectra of 
the magnetic field and of the tension force. 

B3.2 Magnetic-Field Power Spectrum 

There are several formulae that allow one to distill m\ from the 
Stokes correlators. They are all derived by assembling appropri- 
ate linear combinations of the correlators and angle averaging. The 
simplest such formulae are found by noticing that the angle aver- 
ages of the 4th-order parts of E/g, ~L W and of ~Lqq - ILyu vanish [see 
Eq. dB22l l. so 



dtp ReS 



(^-s')mi, 



1 r" 

2n Jo 

_L f 2 " 

2n Jo 

1 r 2 " 

— dtp{-L QQ --L uu ) = SB x B y 



dtpRel, w = 2B x B y m\ 



in 1. 



(B40) 



The standard one-dimensional magnetic-field power spectrum is 
defined with an additional wave-number-space volume factor of 
4nk 2 [Eq. dB6H — the resulting expressions for it are given in 
Eq. {29}. 

It is also possible to construct formulae that formally do not 
require angle averaging at all: using the fact that for certain com- 
binations of the Stokes correlators the 4th-order contributions must 
vanish [assuming isotropy; see Eq. d!6bl. we find 

Re (E ; g sin 2tp - ~L W cos 2tp) 



nit 



"'1 



(b 2 x - B y \ sin 2tp - 2B x B y cos 2tp 
(Sgg - £[/[/) ?,m4tp - 2ReE ef/ cos4</> 
2 ^b] - S^j sin 2tp - 2B x B y cos 2tp 



(B41) 



Since mi is independent of the orientation of the wave vector, these 
expressions can be averaged over the angle tp with arbitrary weight 
functions. Note that in order to compute the shape of the spectrum 
from Eq. dB41t . we do not need to know the magnitude of the mean 
field but we do need its orientation (in the plane perpendicular to 
the line of sight). 

In practice, we expect the formulae given by Eq. dB40t to work 
better than those given by Eq. dB4U because the latter rely on exact 
cancellations that are probably not going to happen with very high 
precision in realistic situations. Even moderate errors in cancelling 
the 4th-order correlators could then easily overwhelm the second- 
order ones: indeed, the terms in Eq. {B39} involving the mean field 

—2 

are small because the mean field was assumed to be weak: B <K 
2 

(b 2 ), so B mi <k Cjj jn „. In Eq. dB40t . the cancellation of the 4th- 
order correlators comes from angle averaging and there is hope that 
mi (k) could be recovered (but see the caveat at the end of Sec. 13.6b . 
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B3.3 Magnitude and Orientation of the Mean Field 



and carrying out the angle averages, we get 



The orientation of the mean field (or, rather, its projection on the 
plane perpendicular to the line of sight) can be easily determined 
from the Stokes parameters themselves: the angle between the 
mean field and the x axis satisfies 



tan 2© 



2B. X B, 



(V) 
<G>" 



(B42) 



Note that tan 2© tells us the orientation but not the sign of the mean 
field. This angle can also be determined from the Stokes correlators 
via Eq. jB40t : 



tan 2© 



dip ReS /(/ J^" d(p (S ee - T, ul} ) 



(B43) 



J^" dp ReS ;e 4 dp ReS ;e 

but the validity of these formulae, unlike that of Eq. dB42t . is pred- 
icated on assuming the statistical isotropy of the fluctuating field. If 
this assumption is well satisfied, the value of tan 20 obtained from 
Eq. dB43t should be independent of k. 

The magnitude of the mean field is a slightly more compli- 
cated quantity to determine. From the total emission intensity [see 
Eq. ©], 



</> = B ± + {K) + (b 2 y ) 



—2 2, 



(B44) 



where B ± = B x + B and the last expression follows from assuming 
the statistical isotropy of the fluctuating field b. Thus, from aver- 
aging /, we can calculate the total energy density of the magnetic 
field but not individually the magnitudes of the mean and fluctuat- 
ing fields. On the other hand, once we know 0, we can find B x Wi 
from Eq. l |B40| > or Eq. dB41t . Let us integrate this quantity over all 
wavenumbers and denote the result by A: 



Jo 



dkAnk l m y (k) = B ± (b ) = A. 



(B45) 



Then {b 2 } = A/B ± and substituting this into Eq. dB44t . we arrive at 
a biquadratic equation for B ± , whose solution is 



2 A 
3<7>' 



(B46) 



B, wi + 4th order, 



(B48) 



S 2 ()t) = 2B' ± mi + 4th order, 
£3 (it) = 4th order, 

3-2 

E4IK) = — B m\ + 4th order, 
4 

where the 4th-order parts are given by Eq. dB23t and real part of 
the Stokes correlators is taken everywhere to eliminate the contri- 
butions from the 3rd-order statistics. 



B3.5 Tension-Force Power Spectrum 

To find the tension-force power spectrum when 6^0, substitute 
Eq. d25t into Eq. i I St and use the solenoidality of the magnetic field 
[Eq. O and Eq. dBTTTl l: 

kj ( c 'ijJ + c y,>) + 

(B49) 

From Eq. dB12l l. we see that the second term vanishes. Using 
Eq. dBlt to express the first term in terms of the magnetic-field 
power spectrum, setting k = £ ± (cos p, sinp, 0) (line-of-sight inte- 
gral) and angle-averaging over p, we get 

G>(jfc) = I k^miik) + 4th order, (B50) 

where the 4th-order part is given by Eq. l |B34| > and is calculated 
from Stokes correlators in the same way as in Sec. IB2.3l Namely, 
using Eq. lB36t and Eq. fB48t . we get 



= '-k L IS, + - ^ 3 + 2Z 4 - \ B 2 ± m 



(B51) 



where Si, . . . , S4 are defined in Eq. dB23t (with real parts taken 
of the Stokes correlators) and B ± mi is calculated via one of the 
formulae given in Sec. IB3.2l Eq. ( 130b follows upon multiplication 
by the wave-number-space volume factor Ank 2 . 

This paper has been typeset from a TgX/ ETgX file prepared by the 
author. 



We have chosen the "-" root because we are assuming B <sc (b 2 ) 

2 

(weak mean field). While we do not really need to know B ± to 
disentangle the second- and 4th-order statistics in Eq. l |B39l l, we 
can use Eq. dB46t to check that the mean field really is weak: 



B ± « (I). 



(B47) 



B3.4 Fourth-Order Quantities 

Now that we know the mean field and mi(k), we can use this infor- 
mation in Eq. flB39t to isolate the 4th-order statistics in the Stokes 
correlators and then proceed to calculating all observable 4th-order 
quantities in the same way it was done in Sec. IB2.2I In general, 
this involves subtracting from the (real part of) the Stokes corre- 
lators the terms that contain m\ [Eq. dB39H so that only the 4th- 
order contributions [Eq. JB22H remain. Doing this requires know- 
ing B ± nii(k) (see Sec. IB3.2l l and the orientation of the mean field 
[Eq. dB42b l. Subtracting the second-order contributions from the 
the averaged Stokes correlators introduced in Eq. dB23t is a partic- 
ularly simple operation: substituting from Eq. dB39t into Eq. dB23t 
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